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ABSTRACT 

We present an observation of the Rossiter-McLaughlin effect for the planetary system 
WASP-3. Radial velocity measurements were made during transit using the SOPHIE 
spectrograph at the 1.93m telescope at Haute-Provence Observatory. The shape of the 
effect shows that the sky-projected angle between the stellar rotation axis and planetary 
orbital axis (A) is small and consistent with zero within 2cr; A = 15t£° deg. WASP-3b 
joins the ~ two-thirds of planets with measured spin-orbit angles that are well aligned 
and are thought to have undergone a dynamically-gentle migration process such as 
planet-disc interactions. We find a systematic effect which leads to an anomalously high 
determination of the projected stellar rotational velocity (vsini = 19.6^ i km s _1 ) 
compared to the value found from spectroscopic line broadening (vsini = 13.4 ± 1.5 
km s _1 ). This is thought to be caused by a discrepancy in the assu mptions made in the 
extraction and modelling of the data. Using a model developed bv lHirano et alJ J2009) 
designed to address this issue, we find v sin i to be consistent with the value obtained 
from spectroscopic broadening measurements (v sin i = 15.7j^3 km s _1 ). 

Key words: stars: planetary systems - stars: individual: WASP-3 - techniques: radial 
velocities - techniques: photometric 



The alignment angle provides an insight into the history of 
these planetary systems. Hot Jupiters are thought to have formed 
far out in the system and have migrated inwards. Migration via 
planet-disk interactions are believed to be relati vely non- violent 
and would not per turb the orbital obliquity dLinet all [l99rl 
iMurrav et alJ 1 1998h . It is thought that the timescale for spin- 
orbit alignment via tidal interactions is of a simi l ar order to the 
lifetim e of the system l lBarker & Qgilviel 120091 : llxvrard et al.l 
2009) and therefore that the primordial alignment of the system 
is preserved. Assuming that a planet was formed in an aligned 
proto-planetary disc, it too is expected to be well alignedQ 
However, more dynamically violent processes such as planet- 
planet scattering and Kozai oscillations from interactions with 
a stellar companion are thought to significantly effect the or- 
bital obliquity and could p roduce highly misaligned an d even 
retrograde orbit (e.g jRasio & FordI 199rj : IWu & Murray||2003l ; 



1 INTRODUCTION 

Bright transiting planets hold the key to determining proper- 
ties of exoplanet systems which are otherwise inaccessible. One 
such example is the projected obliquity of the planet's orbit 
relative to the rotation axis of the star (A). This can be mea- 
sured from spectroscopic observations during transit, which re- 
veal a radial velocity deviation from the standard Keplerian or- 
bital motion known as the Rossiter-McLaughlin (RM) effect 
teossiter|[l924l : lMcLaughlinlll924h . The effect is caused by the 
planet passing over and blocking a portion of the rotating stellar 
surface. This introduces an asymmetry in the spectral line pro- 
file and a small apparent shift in the position of the line. The re- 
sulting shape of the RM effect traces the trajectory of the planet 
across the stellar disc and, once modelled, allows the projected 
spin-orbit alignment to be found. 



* Based on observations collected with the SOPHIE spectrograph 
on the 1.93 m telescope at Observatoire de Haute-Provence (CNRS), 
France, by the SOPHIE Consortium (program 08B.PNP.SIMP). 
f Email: esimpson05@qub.ac.uk 



1 Several authors have noted that this may not necessarily be the case 
in all systems, es pecially those which have undergone int eractions with 
other stars, (e.g.. lGolimowski et alj200dlBate et alj2009h . 
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iNagasawa et al.l l2008). Thus the distribution of alignment an- 
gles allows us to constrain the formation and evolution mecha- 
nisms of planetary systems. 

The form of the RM effect depends on A as well as the pro- 
jected equatorial rotation velocity of the star (v sin i) and sev- 
eral parameters which can be measured from a transit light curve 
(see Section l2~2t , For the WASP-3 system, the amplitude of the 
effect can be estimated as ~ 0.7 v sini 

86 m s -1 , where R p and 7?, are the planetary and stellar radii, 
b is the impact parameter and 0.7 is a scaling factor which ap- 
proximates the effects due to limb darkening. WASP-3 is the the 
fastest rotator for which the RM e ffect has not been p reviously 
reported (v sin i = 13.4 km s -1 . IPollacco et al . 2008). It is also 
a bright target (V^nag = 10.6), and therefore an ideal candidate 
for measuring the RM effect. WASP-3 is a main-sequence star 
of spectral type F7-8V. It hosts a hot lupiter o f mass l.SMj 
and radius 1.3 Rj in a c lose orbit of 1.8 days IPollacco et alJ 
120081 : iGibson et al.ll200l) . Due to the close proximity and rel- 
ative temperature of its host star, WASP-3b is one of the most 
highly irradiated exoplanets known. 

It has b een shown tha t the formulae used to describe the 
RM effect dOhta et alj20 05) appear to overestimate v sin i com- 
pared to the value found from sp ectral line broadening, and the 
devia tion scales with v sini (e.g. , Winn et al ]2005l : lTriaud et al] 
2009). This has led lHirano et al.l d2009l) to re-evaluate the equa- 
tions assuming that the distortion of the radial velocity shift is 
being measured with a cross-correlation technique (as used in 
this work). We apply this method to the RM effect of WASP-3 
an d compare i t to th e results obtained using the model derived 
bv lOhta etaf] ( l2005h . 

In this paper we report on time-series spectroscopic obser- 
vations of WASP-3 obtained during transit. Section[2]describes 
the observations and analysis procedure performed to measure 
the spin-orbit angle of the system. The results of the derived 
system parameters are presented in Section[3] and discussed in 
Section[4] 



2 OBSERVATIONS AND METHOD 

2.1 Observations and data reduction 

A transit of WASP-3b was observed with the SOPHIE spectro- 
graph at the 1.93m telescope at Haute-Provence Observatory on 
the night of 2008 September 30. SOPHIE is a cross-dispersed, 
environmentally-stabilised echelle spectrograph (wavelength 
range 3872.4-6943.5 A) designed for high -precision radial ve- 
locity measurements dBouchy et al.l 120091) . The spectrograph 
was used in high efficiency mode (resolution 7? = 40000) and 
the CCD in fast read-out mode. Two 3 arc-second diameter op- 
tical fibres were used, the first centred on the target and the sec- 
ond on the sky to simultaneously measure its background in case 
of contamination from scattered moonlight. The moon illumina- 
tion was 3% and at a distance of > 86° on the night of transit so 
this did not significantly affect the radial velocity determination. 

We acquired 26 spectra of WASP-3 on the night of transit, 
covering the full transit (137 minutes) and a period of duration 
130 minutes out of transit. The exposure times range from 5 to 
3 1 minutes in order to reach a constant signal-to-noise ratio of 
35 at 550 nm, reflecting the variable throughputs obtained due 
to changing atmospheric conditions. The sky was clear when the 
observations began, and the airmass ranged from 1.1 to 2.7 dur- 



ing the sequence with increasing cloud during the later hours. 
The exposures remained sufficiently short to provide a good 
time sampling with 17 measurements taken during transit. In 
addition, we used 7 out-of-transit observations at various orbital 
phases o btained during t he discovery of the planet in 2007 to fit 
the orbit ( IPollacco et"al]|2008l) . 

The spectra were extr acted and radial velocities measured 
using the SOPHIE pipeline jPerruchot et al.l200 8; Bouchv et alJ 
2009). Radial velocities were computed from a weighted cross 
correlation of each spectr um with a numerical m ask of G2 spec- 
tral ty pe, as described by Baranne et al] d 19961) and lPepe et alJ 
1120021) , A Gaussian was fitted to the cross-correlation func- 
tions to obtain the radial velocity shift, FWHM, and contrast 
with respect to the continuum. The uncertainty in the radial 
velocity was compu ted using the empir i cal re lation given in 
iBouchv et ail d2005l) and ICameron etaT] d2007l) . Uncertainties 
were typically ~ 38 m s _1 during the RM sequence and ~ 23 
m s _1 at other orbital phases. For the observations taken on the 
night of transit, the spectrum below 4367 A was removed from 
the cross correlation in an effort to reduce the systematic errors 
caused by the airmass effect. The wavelength stability over the 
duration of the observations was ~ 3 m s -1 . The 7 orbital points 
reported in Polla cco et al.l d2008l) were re-reduced using an up- 
dated version of the pipeline. Table[T]summarises the measured 
radial velocities. 



2.2 Determination of system parameters 

The RM effect and orbit were fitted simultaneously using all the 
available spectroscopic data. A Keplerian model w as used for 
the or bit, and the analytical approach described in lOhta et alJ 
d2005h (hereafter OTS) for the RM effect. The OTS equations 
were modified to make them dependent on R p /R* and a/R, 
rather than R p , R* and a, to reflect the parameters derived 
from photometry, and reduce the number of free parameters. 
The model comprises 13 parameters: orbital period, P; mid- 
transit time, To ; planetary to stellar radius ratio, R p /R t ; scaled 
semi-major axis a/i?»; orbital inclination, i; orbital eccentric- 
ity, e; longitude of periastron, cj; radial velocity semi-amplitude, 
K; systemic velocity of orbital dataset, 71; systemic velocity 
of transit dataset, 72; sky projected angle between the stel- 
lar rotation axis and orbital angular momentum vector, A; pro- 
jected stellar rotational velocity, v sin i; and stellar linear limb- 
darkening coefficient, u. 

The six observations taken immediately before the transit 
ingress indicate that the systemic velocity of the transit data- 
set is ~ 25 m s _1 lower than that measured from the orbital 
velocities. We therefore allow the systemic velocity of the orbit 
and transit data to differ (71 and 72). All of the out of transit 
points (from 2007 and 2008) were first fitted and no significant 
deviation from a circular orbit was found (e = 0.07±0.08). The 
eccentricity was then fixed to zero for the rest of the analysis. 

We used a linear limb darkening law, as the quadratic law 
alters the model by less than ~ 3ms -1 and so does not seem 
justified by the precision of the RM data, especially given that 
for b ~ 0.5 as in this case, the alignment angle is primar- 
ily determined by the time when the radial velocity anomaly 
vanishes, rather than the exact form of the model. The linear 
limb-da rkening coeffi cient u = 0.69 was chosen from the ta- 
bles of lClareJ d2004h (ATLAS models) for the g' filter . We 
used the value corresponding to T e g — 6500 K, logg = 4.5, 



© 0000 RAS, MNRAS 000, 000-000 



The spin-orbit alignment of WASPSb 3 



Table 1. Radial velocities and 1<t error bars of WASP-3 measured with 
SOPHIE during and outside transit. The orbital points are updated from 
IPollacco etaflfeoOgl) using a new version of the reduction pipeline. 



BJD 


RV 


Error 


-2 400 000 


(km s _1 ) 


(km s _1 ) 


Planetary transit: 






54740.29509 


-5.42487 


0.03856 


54740.29909 


-5.38158 


0.03628 


54740.30328 


-5.36639 


0.03751 


54740.30771 


-5.43318 


0.03734 


54740.31312 


-5.45738 


0.03756 


54740.31830 


-5.40972 


0.03828 


54740.32351 


-5.42139 


0.03816 


54740.32907 


-5.33800 


0.03782 


54740.33524 


-5.35707 


0.03764 


54740.34152 


-5.38229 


0.03741 


54740.34699 


-5.38596 


0.03626 


54740.35300 


-5.36107 


0.03533 


54740.35891 


-5.36891 


0.03508 


54740.36506 


-5.41906 


0.03437 


54740.37170 


-5.49966 


0.03580 


J4 /4U. j llli 


-j.jjjUj 


U.UJooo 


54740.38426 


-5.55488 


0.03571 


54740.39228 


-5.58417 


0.03825 


54740.39993 


-5.64022 


0.03606 


54740.40758 


-5.58523 


0.03519 


54740.41466 


-5.69115 


0.03415 


54740.42259 


-5.64015 


0.03543 


54740.43066 


-5.53109 


0.03945 


54740.44171 


-5.59288 


0.03844 


54740.45894 


-5.62403 


0.04132 


54740.48062 


-5.52085 


0.05355 


Other orbital phases: 






54286.52255 


-5.72553 


0.02137 


54287.45639 


-5.18430 


0.03541 


54289.36622 


-5.18792 


0.03047 


54340.32511 


-5.64313 


0.01280 


54341.39898 


-5.40394 


0.01914 


54342.31981 


-5.54401 


0.02801 


54343.48259 


-5.60354 


0.01283 



[M/H] = 0.0 a nd C = 2 km s x (a typical value for this 
spectral type, see iGrav et alj|200lh which most closely repre- 
sented the pro perties of WASP-3 (T e ff = 6400 K, log g = 4.25, 
[M/H] = O.O. IPollacco et al.ll2008h . As a test, we set u as free 
parameter and found that there was no significant effect on A or 
v sin i. 

We placed no constraint on v sin i as there is a known sys- 
tematic effect in its determ i nation from RM measurements (see 
Section [3~2l l. iTriaud et alj (2009) suggest that there may be a 
similar systematic effect influencing R p /R, if solely spectro- 
scopic observations are usedQ We therefore fixed the photo- 
metric parameters (P, To, R p /R„, a/R, and i) to those found 
through light curve modelling. An MC MC fit to a full an d partial 
transit of WASP-3 was presented in lGibson etahl d2008l) . In this 
analysis, the stellar radius was calculated based on the scaling 

1/3 

relation R, oc M* with a prior on M* , leading to an underes- 

2 This is because the amplitude of the RM effect is primarily deter- 
mined by Rp/R* and vsini, so the observed increased amplitude 
could lead to a false inflation of these values. 



Table 2. Adopted system parameters and uncertainties requir ed for the 
RM effect, and other photometric parameters updated from Gi bson et alj 
j2008h . 



Parameter (units) 


Symbol 


Value 


Period (days) 


P 


1.846835 


± 0.000002 


Transit epoch (HJD) 


To 


2454605. 55915± 0.00023 


Orbital inclination (deg) 


i 


84.93 


+ 1.32 
-0.80 


Planet/Star radius ratio 


Rp/R* 


0.1013 


+0.0014 
-0.0013 


Scaled semi-major axis 


a/R, 


5.173 


+0.246 
-0.162 


Eccentricity * 


e 





(adopted) 


Limb darkening ' 


u 


0.69 




Impact parameter 


b 


0.466 


+0.047 
-0.111 


Transit duration (hours) 


T d 


2.737 


+0.024 
-0.022 


Stellar radius (Rq ) 


R* 


1.31 


+0.05 
-0.07 


Planet radius (Rj) 


Rp 


1.29 


+0.05 
-0.07 


* See Section[2] 








t from Claret (2004) 









timation in the error on the orbital inclination. The light curves 
were re-fitted with the stell ar radius as a free p arameter using 
the procedure described in iGibson et al] (2009). Five separate 
chains each with two-hundred thousand points were computed. 
To test that the chains had all converge d to the same parameter 
space, the Gelman & Rubin statistic dGelman & Rubinlll992l) 
was calculated for each of the free parameters. It was found to 
be less than 1% from unity for all parameters, a good sign of 
mixing and convergence. The results are shown in Table[2] 

Thus the model contains 5 free parameters (A, v sin i, K, 
71 and 72). The y 2 statistic was calculated: 

* - L [ j t J « 

i 

where Vj lC> bs and i>i iCa i c are the ith observed and calculated ra- 
dial velocities and a% the corresponding observational error. The 
optimal parameters were obtained by minim ising the y statis- 
tic using the IDL-based MPFIT function (Markwardti l2009h : 
a least-squares minimisation technique using the Levenberg- 
Marquardt algorithm. The starting values for the parameters 
were those derived in Pollacco et al. (2008) as well as A = 
and 72 = 71 . 

The spectroscopic errors were rescaled so that the best 
fitting model has a reduced y 2 (Xred) °f 1> which required 
the initial spectroscopic errors to be multiplied by the factor 
yj Xr e d ~ 0.905. This suggests that the errors were slightly 
overestimated however the rescaling did not significantly affect 
the best-fit parameters or errors. 

The la best-fit parameter uncertainties were calculated us- 
ing a Monte-Carlo method. We created 10° synthetic data sets 
by adding a la Gaussian random variable to the data points. 
The photometric parameters were also varied for each realisa- 
tion by taking a parameter set from the distributions from the 
light curve fitting. This allowed correlations between parame- 
ters such as a/R, and i to be taken in to account. The starting 
values for the free parameters were also perturbed by a 3a Gaus- 
sian random variable to further explore the parameter space. The 
free parameters were re-optimised for each simulated data-set 
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Figure 1. Upper panel: Radial velocities of WASP-3 minus the sys- 
temic velocity, overp lotted with the best-fit model. Open points are from 
iPollacco et at] J2008h . Filled points show the RM effect of WASP-3 on 
2008 September 30. Lower panel: residuals from the best-fit curve. 
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Figure 2. Upper panel: RM effect of WASP-3 minus the Keplerian mo- 
tion, overplotted with the best-fit model. The solid and dotted lines show 
the offset between the radial velocity and photometric zero-points, sug- 
gesting that the system is not completely aligned. The error bar on 
the solid line shows the lrr uncertainty in the position of the radial 
velocity zero-point. The uncertainty on the photometric zero-point is 
smaller than the width of the line. Lower panel: residuals from the best- 
fit model. 



to obtain the distribution of the best-fit parameter values. The 
distributions were not assumed to be Gaussian and the la lim- 
its were found where the distribution enclosed ±34.1% of the 
values away from the median. In order to verify the errors, we 
also ran a global fit to both the photometric a nd radial veloci- 
ties u sing an MCMC algorithm, as described in dAnderson et al] 
l2009h . 



3 RESULTS 

3.1 System parameters 

The derived parameters and errors for the WASP-3 system are 
shown in Table [3] for both the least-squares and MCMC meth- 
ods. The results from the global MCMC analysis are consistent 
with those found from the least squares fitting. The re-reduced 
radi al velocities yield or bital parameters which are consistent 
with lPollacco et al] ( feOOSt . 

Figure[T]shows all the WASP-3 radial velocities minus the 
systemic velocity; the best-fit model; and the residuals to the 
fit. Figure [2] shows a zoom-in on the observations taken on the 
night of transit minus the Keplerian orbit; the best fit model; 
and residuals. It can be seen that there is a ~ 10 minute offset 
between the radial velocity and photometric zero points (solid 
and dotted lines). The elapsed time between the transit photom- 
etry and the RM observations is approximately 73 orbits and the 
0. 173s uncertainty in P and 20s uncertainty in To causes an un- 
certainty of only 33 s which cannot explain the offset. By fitting 
the data, we obtain A = 15i 9 ° deg, indicating that the timing 
offset is caused by a small misalignment between the planetary 
orbit and stellar spin axis. The system is generally well aligned 
and the misalignment is not highly significant, with A consistent 
with zero to within 2a. 

The difference between the offset velocities of the 7 orbital 
data points and those taken on the night of transit is significant 
to 2 a; ~ 28 ± 14 ms _1 , so we investigated modelling all the 
data together without including a relative shift. We found that 
the newly derived parameters were consistent with the previ- 
ous values to within the error bars. The value of A increased 
to 23ly deg because A is highly correlated with 7. The new 
fit yields a 7 velocity which is higher than that found for the 
transit dataset alone, so the RM effect appears to cross the zero 
velocity point earlier, and produces a larger misalignment angle. 
There are several factors which could cause an offset between 
the datasets. These include a long-term trend caused by an addi- 
tional body, stellar activity, instrumental drift or a combination 
of these In addition, the data taken on the night of transit un- 
derwent a slightly different reduction procedure which may well 
cause a shift in the radial velocities between nights, (see Section 
[2} so we prefer to fit the two datasets separately. 

3.2 Stellar rotational velocity and systematics 

The observed amplitude of the RM effect (~ 100 m s _1 ) 
is significantly larger than expected (~ 86 m s" 1 , see Sec- 
tion [TJ. Figure [3] shows the predicted RM effect with v sin i 
equal to that found from modellin g the spectral line broadening 
(v sini = 13.4 ± 1.5 km s~ ^ IPollaccoetaLboOSh . By fitting 
the RM effect, we find the projected rotational velocity to be 
overestimated by ~ 50%; wsini = 19.6tH Ions' 1 . Several 
authors have noted a similar discrepancy, which is particularly 
evident in faster rotating stars. IWinn et al.l J2005h suggested that 
systematics are present in the determination of the radial veloc- 
ities because the asymmetric stellar line profile is being fitted 
with a symmetric Gaussian. They attempted to correct this by 
adding a quadratic term to the RM model which significantly 
improved the estimated v sin i. 

3 The shift is larger than the typical value for astrophysical jitter for 
this spectral type so this is unlikely to be the sole cause. 
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Table 3. Best-fit values and uncertainties using Ohta et al. 1 2005) (OTS) (both Least-squares and MCMC fitting) and lHirano et al] f2009) (H09) models. 



Parameter (units) 


Symbol 


OTS 


OTS (MCMC) 


H09 


Projected alignment angle (deg) 


A 


IB ij? 


14 ±9 


13 1? 


Projected stellar rotation velocity (km s -1 ) 


v sin i 


19.6 1 2 ,;? 


19.6 ±2.2 


15-7 


RV semi-amplitude (km s — 1 ) 


K 


0.274 ± 0.011 


0.273 ± 0.013 


0.276 ± 0.011 


Systemic velocity of orbit dataset (km s — 1 ) 


71 


-5.458 ±0.008 


-5.458 ± 0.005 


-5.458 ± 0.007 


Systemic velocity of transit dataset (km s ) 


72 


-5.483 ±0.011 


-5.486 ± 0.011 


-5.487 ± 0.009 


Chi-squared (spectroscopic data) 


x 2 


27.9 


27.4 


27.2 



iTriaud et alj d2009l) showed that in HD 189733, the struc- 
ture of the residuals observed when subtracting the RM model 
from the data was consistent with those found when subtract- 
ing the RM model from simulated data undergoing the same 
data reduction procedure. This indicates that the cause of the 
discrepancy is indeed in the extraction of the radial velocities. 
The residuals of the WASP-3 RM effect show a very similar 
structure d uring the first three quarters of the transit to those in 
Figure 4 of Tria ud et al.1 d2009l) (the last section is dominated by 
systematics caused by airmass and weather effects). The ampli- 
tude of the residuals is of the order of the error bars so it is not 
clear whether the structure is detected significantly. 

The apparent radial velocity shift of the stellar lines pre- 
dicted by the OTS model is produced by the distortion of the 
line profiles by the planet blocking of a velocity component of 
the star's spectrum. Assuming that the resulting asymmetry of 
the stellar lines is not well resolved, a symmetric function fit- 
ted to the line profile measures a shift in the mean position. For 
faster rotating stars, the asymmetry in the line profile is better 
resolved because of the broader width, and causes a larger ap- 
pare nt shift than expecte d, as seen here. 

Hirano et aD {2009) (hereafter H09) have addressed this 
problem by modifying the OTS equations under the assumption 
that the RV shift is measured by cross-correlation with a Gaus- 
sian function. They find that the RM model is also dependent on 
the intrinsic width of the individual line profile. 

Hence, instead of modelling the RM velocity anomaly us- 
ing the OTS equation; 



Sv = fv p /(l - /) 



(2) 



where / is flux blocked by planet ~ (R p /R*) 2 and v p is sub- 
planet velocity (ie the velocity component of the rotating stellar 
surface blocked by the planet ~ x p v sin i /R*, where x p is the 
x co-ordinate of the position of the planet on the stellar surface, 
see Fig 5 of OTS), we now use the H09 model (Equations 41 
and 42): 



Sv = fv p {p - qv p ) = fv p (1.72 - 0.00546 v^) 



where 



1 + 



2/3 2 + a 2 



3/2 



2/3 2 



(3) 



(4) 



and P is the intrinsic line width, which has been calculated by 
Miller et al. (in prep) for WASP-3 as /3 = 11.1 ± 0.1km s" 1 . 
The quantity a = v sin i/a, where a is a scaling factor depend- 
ing on limb darkening parameters (see Equation F6 of H09). 
We calculate a = 1.31 for the fixed limb darkening coefficients 
Mi = 0.69 and U2 = 0. 



We fit this equation to the data and find that A is con- 
sistent with that found using the OTS model, A = 13^y deg 
(A = 15tg° deg; OTS). This is expected because A mainly de- 
pends on the shape of the waveform and the time crossing zero 
point rather than the amplitude. The slight difference is proba- 
bly caused by fit finding a marginally lower value of 72, which 
is correlated with A. We find vsmi to be much closer to the 
line-broadening value vsmi = 15.3 ± 1.4 km s _1 . It is still 
slightly different but is acceptable considering the errors. There 
was no significant difference in \ 2 between the OTS and H09 
fits. Figure [3] shows the data over-plotted with the best-fit OTS 
and H09 models. 

We also investigated whether the discrepancy between the 
two independent m easurements of v sin i could be caused by 
differential rotation. iGaudi & Wind 00071) show that if the sur- 
face speed of a star is dependant on latitude as well as longitude, 
then the form of the RM effect is slightly altered and depends on 
the fractional difference in rotation speed between the pole and 
equator (a). To reproduce the ~6 km s -1 diff erence in v sin i 
seen in WASP-3, a must be greater than 0.45. lReinersl l l2006l) 
showed that no stars have been observed with a > 0.45 so this 
scenario seems unlikely. In addition, the form for the differen- 
tial rotation law would have to be anti-solar, which is seldom 
observed. The bias caused by the fitting of gaussians to asym- 
metric CCFs offers a much more viable explanation. 



4 CONCLUSIONS 

A full transit of WASP-3b was observed with the SOPHIE spec- 
trograph at the 1.93m telescope at Haute-Provence O bservatory. 
We m odelled the RM effect using the equations of lOhta et all 
and found the sky projected spin-orbit alignment angle 
of the system to be A = 15Zg deg. WASP-3b is not strongly 
misaligned and is consistent with 0° to within 2a. This sug- 
gests that WASP-3b underwent a relatively non-violent migra- 
tion process which did not perturb it from the primordial align- 
ment of the proto-planetary disc. The result indicates a marginal 
detection of a small misalignment, and further, higher precision 
observations of this target would establish the reality of this. 

We find a significant discrepancy between the v sin i value 
from the RM model (vsmi = 19.6!?, 1 km s_1 ) an d the mea- 
surement from spec troscopic line broade ning analysis (v sin i = 
13. 4± 1.5 km s~ 1 ; |Pollacco et alj 2008). This phenomenon has 
been noted in several systems, especially for fast rotators. It is 
thought to be the result of systematics introduced when measur- 
ing a line profile asymm etry as a shift in the mean line position. 

iHirano et al.1 {2009) provide a modified formulation to take 



© 0000 RAS, MNRAS 000, 000-000 



6 E. K. Simpson et al. 



-200 



100 


-100 




-0.04 -0.02 0.00 0.02 0.04 0.06 

Orbitol phase 

Figure 3. Upper panel: RM ef fect of WASP-3 w ith various models. 
Solid (red) is the fit using the iHirano et all <2009l) (H09) model with 
v sin t as a fre e parameter (usini = 15.3 km s — 1 ). Dotted (blue) 
is the Iphta et al J 120051) (OTS) model with vsini as a free parameter 
(v sin i = 19.3 km s~ ). Dash dot (blue) is the fit using the OTS model 
but with v sin i held fixed at v s ini = 13.4 km s~ 1 as determined 
through spectroscopic broadening i Pollacco et al. 2008). The amplitude 
of this model is obviously too low. All the models cross the x-axis at the 
same phase, giving consistent values of A. Lower panel: residuals from 
the H09 best-fit model. 



into account the method used to determine radial velocity shifts. 
We applied this model and found the alignment angle to be con- 
sistent with the OTS model, and derived an improved value for 
the rotational velocity v sin i = 15.7t};3 kms" 1 . This is con- 
sistent with the value found from spectroscopic line broadening 
analysis. The RM effect has the potential to provide a much 
more precise determination of v sin i than spectral broadening 
for slow rotators, as it is a geometrical measurement that is not 
subject to uncertainties such as macroturbulence. However, RM 
models must be shown to ensure that the v sin i discrepancy is 
fully resolved. 

Since the beginning of 2009, the number of systems re- 
ported with tilts significa ntly greater than 30 ° has risen dra- 
m atically from one (XQ3, lHebrard et al]|2008l) to six (CoRoT- 
l. Jgont et al.ll2009al ; HAT-P-7, |Wir - 1 1 
Moutou et al] 



Pont etal.1 [2009a; HAT-P-7. Iwinn et alj|2009bl: iNarita et alj 



20091: 



HD 80606, 



2009; Pontetal. 



Gillonll2009l;lwinn et all2009al; W ASP-14. |johnson et alj2009l ; 
WASP-17, lA.nderson et alj|2009l) . The RM effect has now been 
measured for a total of 18 systems and WASP-3 joins the ~ 2/3 
of these which are well align ed. These statis tics are re-enforcing 
the inference by iFabrvckv & Winnl J2009I) that there is a bi- 
modal distribution of spin-orbit angles, with a fraction of sys- 
tems being well aligned while the rest have random mutual ori- 
entations, suggesting that there are multiple migration mecha- 
nisms at work. It is therefore vital to increase the number of 
measured systems, both aligned and misaligned, to place these 
theories on a secure statistical basis. 
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